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Several systems involve spatial arrangements of elements such as molecules or cells, the character- 
ization of which bears important implications to biological and physical investigations. Traditional 
approaches to quantify spatial order and regularity have relied on nearest neighbor distances or the 
number of sides of cells. The current work shows that superior features can be achieved by consid- 
ering angular regularity. Voronoi tessellations are obtained for each basic element and the angular 
regularity is then estimated from the differences between the angles defined by adjacent cells and 
a reference angle. In case this angle is 60 degrees, the measurement quantifies the hexagonality of 
the system. Other reference angles can be considered in order to quantify other types of spatial 
symmetries. The performance of the angular regularity is compared with other measurements in- 
cluding the conformity ratio (based on nearest neighbor distances) and the number of sides of the 
cells, confirming its improved sensitivity and discrimination power. The superior performance of 
the haxagonality measurement is illustrated also with respect to a real application concerning the 
characterization of retinal mosaics. 

PACS numbers: 05.50.+q, 64.60.Cn, 87.80.Vt, 05.65.+b 



I. INTRODUCTION 

Several important properties of biological systems are 
directly related and even determined by the spatial dis- 
tribution of their constituent elements. For instance, the 
distribution of ganglion cells through the mammals retina 
is known to accurately sample the visual field with just 
the right amount of overlap. Similarly, the adjacency of 
cells is important for cell signaling, while the spatial dis- 
tribution of cells actively expressing some specific gene is 
immediately related to tissue and organ formation. Other 
examples where the spatial distribution of the elements 
is critical involve the spread of pathological agents such 
as viruses and bacteria as well as the spatial arrangement 
of diseased cells and lesions. 

Typically, spatial order in biological systems involves 
structured/symmetric arrangements of points such as in 
hexagonal systems. These structures can be understood 
as regular lattices, characterized by fixed spacing and 
angles between the constituent points. Even in the cases 
where the boundaries between the elements are not avail- 
able, they can be obtained from the Voronoi tessellation 
considering the original points as seeds (e.g. 53 )■ The 
ubiquity of polygonal organization in biological systems 
is related to special packing and physical properties al- 
lowed by ordered systems (e.g. For instance, it is 
known that the hexagon is the regular polygon with max- 



imum number of sides which can be used to tile the plane. 
The immediate advantage of such a regular tiling is the 
maximization of the number of neighbors of each element, 
with immediate implications for cell signaling : 8| , resolu- 
tion and isotropy in sensory acquisition (e.g. photorecep- 
tors, retinal neurons and omatidea), structural properties 
(e.g. choral, cell membrane, striate muscle fibers, and 
other structures), to name but a few cases. Interesting 
combinations of pentagons and hexagons are also often 
found in three-dimensional biological structures such as 
pollen, viruses and radiolaria structures. 

Because spatial regularity plays such a key role in 
defining specific structural and functional properties of 
biological systems at micro, meso and macroscopic spa- 
tial scales, it becomes of paramount importance to quan- 
tify in an objective and accurate way the degree of spa- 
tial organization of such systems. Among the previous 
works aimed at automated quantification of spatial distri- 
butions, the most frequently adopted approach involves 
counting the number of sides of each cell 0, ^| . An- 
other traditional approach which can be used to quan- 
tify polygonality involves the use of the distance be- 
tween nearest neighboring points as well as derived mea- 
surement such as the conformity ratio (e.g. |3. ITEl IT^|). 
A comparison of methods for characterization of spa- 
tial properties of retinal mosaics has been presented by 
Cook !3|, with emphasis on nearest neighbor distances. 
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The present article starts by identifying the main prop- 
erties required from a good measurement of spatial order 
and follows by describing each of the considered measure- 
ments, whose properties are then evaluated for the char- 
acterization of global and local spatial order. The analy- 
sis of the potential of the measurements to characterized 
global properties takes into account simulated (hexago- 
nal lattices with progressive perturbations) as well as real 
data (retinal mosaics of photoreceptors) . The possibility 
to use the measurements to quantify local order around 
each point is also addressed with respect to the identifi- 
cation of regions with distinct spatial order and the anal- 
ysis of a system involving progressive (radial) variation 
of spatial order. 



II. REQUIREMENTS 

Before considering possible indicators of spatial order 
and evaluating their performance, it is important to iden- 
tify the main features which could be expected from a 
good such measurement. Such features are discussed in 
the following. 

Let us consider that the property P to be measured is 
a function of a given parameter s IS], so that P = P(s), 
and that the act of measuring P involves mapping it into 
the quantity Q(P(s)) — Q(s). For instance, s can be 
the intensity of spatial perturbation added to a perfectly 
hexagonal lattice (see Section HV A|) . and Q is the mea- 
surement of spatial order. Note that typically we do not 
have access to the values of P, otherwise there would be 
no need to consider its measurement. Figure^illustrates 
mappings from the property P into two possible mea- 
surements Q(s). Although this example assumes that 
P is upper and lower bound, some situations may im- 
ply unlimited values of P m in or P max . Situations such 
as that depicted in Figure ^ can be always normalized, 
by shifting the function along the x- and/or y-axes, such 

that P m in — Qmin — and P m ax — Qmax = lj defining 
a linear relationship. The mapping in (b) can also be 
linearized by using an additional transformation of the 
curve itself. 

Sensitivity /Linearity: One important feature of any 
measurement is its ability to respond to variations of 
the measured quantity. Because measurements can be 
normalized, sensitivity should be considered in relative 
terms along the interval of measurement of interest. For 
instance, the mapping in Figure ^a) is more sensitive 
than that in (b) for small values of P(s), but becomes 
less sensitive for large values of P(s). The sensitivity of 
the measurement in Figure^a) is the same for any value 
of P(s) . Observe that the increase of sensitivity along one 
of the regions of the mapping can only be achieved at the 
expense of loss of sensitivity at another region. Unless 
some specific interest of achieving enhanced sensitivity 
at some specific range of values is to be considered, mea- 
surements should ideally be linear. In the specific case of 
spatial order quantification, one is frequently interested 
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FIG. 1: Two possible measurements mapping the property P 
into Q. 



in the cases of small deviations from perfect order. In 
such cases, it is interesting to obtain greater sensitivity 
for higher values of spatial order. 

Generality, Uniformity and Invariance: It is usually 
interesting that the measurement of P be monotonic (in- 
creasing or decreasing), well-defined and with uniform 
properties along the whole interval of measurements. 
This implies absence of singularities (e.g. tendency to 
go to infinity for some values of P) as well as uniform 
sensitivity. The linear mapping meets all such require- 
ments. In cases such as the quantification of spatial or- 
der of geometrical systems, it is also important that the 
measurements be invariant to translation, rotation and 
scaling. 

Discriminative Power: Although Figure ^ considered 
a deterministic mapping from P into Q, it often happens 
that measurements taken over different realizations of a 
system for the same parameter s yield different values. 
This is often the case with spatial order. For instance, 
the spatial order of the photoreceptors in the retina of 
two similar animals (i.e. same age, gender, etc.) from 
the same species will unavoidably be characterized by 
distinct, even though similar, values. In turn, such a 
stochastic variation of Q implies that in certain situa- 
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tions two realizations of a system, respective to different 
parameters s, will yield the same measurement value. 
In other words, the stochastic variability of Q will tend 
to limit its discriminative power. Given two instances 
of s, i.e. si and S2, the discriminative power will vary 
inversely with the overlap between the statistical den- 
sity functions characterizing -P(si) and P(s2). Such an 
overlap can be estimated by considering the average and 
standard deviation of those two densities, which is the 
procedure adopted in the present article. 

Adimensionality: While this is not an imperative prop- 
erty, it is interesting that the measurement in question 
be adimensional, allowing it to be expressed in percent- 
age values. In addition, adimensionality often facilitates 
multivariate statistical analysis and testing (e.g. |9lll0j). 



III. METHODS 



Let the image containing the biological system under 
analysis be represented in an N{ x Nj image A. Each 
element of the system is identified by an integer label i, 
while its spatial position is indicated in terms of the coor- 
dinates (x(i),y(i)) of some reference point or seed, often 
corresponding to the centroid of the elements (see Fig- 
urcEt). The nearest neighbor distance between a point i 
and a set of points S is defined as the smallest distance 
between i and each of the points in S. The Voronoi tes- 
sellation of the system of reference points can be obtained 
by assigning to each pixel A(i,j) the label of the nearest 
reference point. One of the nice features of Voronoi tes- 
sellations is that they can be used to establish boundaries 
and adjacencies between the original reference points. 
More specifically, given a reference point k, each of its 
Voronoi neighbors will correspond to a reference point i 
such that the Voronoi cells of k and i share a common 
side. In order to avoid border effects, the Voronoi cells 
which are adjacent to the image background (i.e. the 
border cells) are henceforth identified and excluded from 
the calculations. 

Number of Voronoi cells: An interesting measure- 
ment which can be immediately obtained from Voronoi 
tessellations corresponds to the number of sides of each 
Voronoi cell, a measurement which has been called hexag- 
onality (e.g. [rHH^ I. Henceforth nk represents the num- 
ber of sides of the cell associated to the reference point 
k. An hexagonal system, for example, will produce all 
cells with 6 sides, except for the cells at the border of the 
system. 

Angular Regularity: For each Voronoi cell k, start- 
ing at an arbitrary neighbor i, determine the successive 
Nk neighbors as one turns around the cell in a clockwise 
fashion and calculates the angles oti defined by the ad- 
jacencies (see Figure Et)- The angular polygonality of a 
given point k can then be expressed in terms of the sum 



of angle differences for each point, i.e. 



Ei 

i=l 



(1) 



where [3 is a specific angle of interest. In case one is in- 
terested to quantify the angular hexagonality of a point 
system, j3 should be taken as 60 degrees. Similarly, in 
case one wants to check for orthogonal arrangement, (3 
should be 90 degrees. Other values of (3 varying in the 
interval < (3 < 180° can also be used in order to check 
for other types of angular distributions. However, only 
60, 90 and 120 degrees allow full regular tiling of the 
plan. Figure |3 shows the polygonalities of a single point 
surrounded by N points (uniform angle spacing). Note 
that the sum of angle differences varies from to infinite 
(i.e. a point surround by an infinite number of points). 
Although it would be possible to consider the average of 
angle differences (i.e. divide the sum of angle differences 
in Equation^by the number of neighbors of each point), 
such an alternative would lead to similar measurements 
as the number of neighboring cells increase, therefore re- 
ducing the discriminative power. 

An alternative measurement of spatial order, yield- 
ing the adimensional measurement henceforth called the 
polygonality index, can be defined as a function of the 
sum of difference angle as follows 



A a (fc) 



1 



(2) 



This measurement, which was first considered for 
hexagonality characterization in jj], is bound between 
(lack of spatial order) to 1 (perfect polygonality). 

Conformity Ratios: Once all nearest neighbor dis- 
tances have been calculated, for instance by using the 
simple and effective algorithm described in [3(, the con- 
formity ratio of distances can be defined |l6j as the ratio 
between the mean and standard deviation of all near- 
est neighbor distances. More organized systems should 
therefore imply higher conformity ratios. 

Coefficient of Variation: The exact inverse of the 
conformity ratio, the coefficient of variation is also con- 
sidered in this work because of its traditional use in 
statistics. 

Note that the conformity ratio of the nearest distances 
present a singularity (tends to infinity) when the short- 
est distances for every point in the system is equal, as 
with regular lattices. The coefficient of variation does 
not present singularities, but will produce null values for 
any regular lattice (null standard deviation of shortest 
distance will imply coefficient of variation equal to zero). 

Although it would be possible to define conformity ra- 
tio and coefficient of variation of angles differences, such 
measurements would be completely unable to cope with 
situations such as that considered in Figure^ i.e. involv- 
ing N equal (or similar) angles around a point. In these 
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(a) 



(b) (c) 

FIG. 2: The geometrical construction and symbols used in 
the definition of the angular hexagonality are presented in 
(a). The reference points are shown as black dots, and the 
neighbors of seed k are enumerated in clockwise fashion from 
1 to Nk- Although rather distinct, the two point distributions 
in (b) and (c) yield identical statistics of shortest distances, 
and hence the same conformity ratio. 



cases, the standard deviation of the angle values is null, 
implying conformity ratio equal to zero and coefficient of 
variation equal to infinity Such a property also implies 
lack of sensitivity for systems characterized by a high 
level of spatial order. Contrariwise, the sum of angle dif- 
ferences or hexagonality index will produce distinct val- 
ues for any number of equally spaced neighboring points. 
For such reasons, and also because of poor performance 
in preliminary experimental investigations conducted by 
the authors, the coefficient of variation and conformity 
ratio of angle differences will not be considered further 
in this article. 

Although all the above measurements are invariant to 
translation, rotation and scaling of the system under 
analysis, any measurement derived from the minimum 
distance statistics will fail to distinguish between situ- 
ations such as those illustrated in Figure Efb) and (c), 
which yield identical statistical densities of shortest dis- 
tances and therefore equal conformity ratios. This im- 
portant limitation of the shortest distance measurement 
stems from the fact that it only takes into account the 
most immediate neighborhood around each point, over- 
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FIG. 3: The poligonalities for /3 = 120, 90 and 60 degrees 
obtained for a single point surrounded by N uniformly dis- 
tributed angles (i.e. a — 2n/N). Observe that the sharp 
peaks clearly identify each type of symmetry. For instance, 
the perfect hexagonal arrangement of neighbors for TV = 6 
implies the sharp peak of the poligonality assuming [3 — 60 
degrees at that respective position. 



looking spatial properties at higher spatial scales. There- 
fore, a natural means to try to obtain better discrimina- 
tive power while using distance statistics consists in con- 
sidering successive shortest distances, such as the second, 
third, etc., shortest distances. Such a possibility is also 
explored in this work. 

A severe limitation of the method of counting the num- 
ber of Voronoi cell sides follows from the fact that the 
number of cell sides will not change except for gross per- 
turbations of the particle system. However, this mea- 
surement can be potentially useful when the structures 
under analysis involve large perturbations. 

IV. RESULTS 

This section describes a systematic performance com- 
parison of methods used to quantify the spatial order of 
point distributions, including nearest distances, number 
of Voronoi sides, conformity ratio and coefficient of varia- 
tion of the nearest distance, sum of angle differences and 
hexagonality index. This evaluation considers the quan- 
tification of global and local spatial order in systems of 
points. While the former case addresses the problem of 
assigning a single measurement to the whole system of 
points, the latter situation involves the quantification of 
the spatial order around each point of the system. The 
global evaluation is performed for simulated data (regular 
hexagonal lattices with progressive spatial perturbations) 
and real data concerning agouti (Dasyprocta agout) reti- 
nal photoreceptor mosaics. The obtained results clearly 
corroborate the superiority of the sum of difference angles 
and the hexagonality index, especially regarding their lin- 
earity, sensitivity and discriminative power. The inves- 
tigation of the potential of the considered measurements 
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for quantifying local spatial order involves the identifica- 
tion of regions with different levels of order as well as a 
system where the spatial order varies radialy. 



A. Simulated Data 



In order to investigate the performance of the above 
measurements under a controlled situation, hexagonal 
lattices with progressive spatial perturbation intensi- 
ties were obtained by using mathematic-computational 
means. More especihcally, hexagonal lattice with 621 
points and vertices of length A = 10 pixels were gener- 
ated in a 480 x 480 pixels image and progressively pertur- 
bated with uniformly distributed displacements of mag- 
nitude 6 = 1, 2, . . . , 10. A total of 50 realizations of each 
configuration (i.e. each perturbation intensity) was per- 
formed in order to enhance statistical representativity. 
Figure^lshows two of the considered simulated hexagonal 
lattices, respective Voronoi tessellations and hexagonal- 
ity indices obtained for perturbation intensities 6 = 2 and 
4, respectively. The histograms of the number of Voronoi 
cell sides, coefficient of variation of shortest distances, 
conformity ratio of shortest distances, shortest distance 
values, average angle difference (in degrees), and hexago- 
nality indices for each of the perturbation intensities are 
shown in Figure El 

It is clear from these results that the number of Voronoi 
cell sides, shown in Fig. Efa), is largely invariant to the 
perturbation intensity. In the cases where the latter is 
large enough to induce changes, it acts mainly by in- 
creasing the dispersion of the measurements while keep- 
ing an almost constant average value. The coefficient 
of variation (b), first shortest distance (d) and sum of 
angle differences (e) all resulted in an almost linear map- 
ping between the perturbation intensity and the mea- 
surement. All these cases, and especially the coefficient 
of variation, are also characterized by an increase of dis- 
persion for larger perturbations. The conformity ratio re- 
sulted the most non-linear relationship, providing greater 
sensitivity for smaller perturbations. However, the mea- 
surement dispersions are also substantially larger in such 
cases, which limits the discriminative power of this mea- 
surement. The hexagonality index, shown in Fig. Elf), 
accounts for an almost linear mapping, with a slight in- 
creased sensitivity for small perturbations. In addition to 
these interesting features, this measurement is the only 
one characterized by an almost constant standard devi- 
ation which, combined with the almost linear mapping 
characterizing this measurement, implies a constant dis- 
criminative power along all the considered perturbation 
intensities. The linearity, sensitivity and discriminative 
power of the measurements of spatial order are summa- 
rized in Table |U It follows from such results that the 
hexagonality index is particularly suitable in cases of 
higher spatial order, while the sum of angles provides 
an interesting alternative for analysis of less ordered sys- 



tems. 

The properties of the number of sides, shortest distance 
and hexagonality index are confirmed by analysis of the 
histograms in Figure El which show the distribution of 
such measurements for several perturbation intensities. 



B. Biological Data 



In order to further evaluate the performance of the 
conformity ratio and sum of angle differences, these two 
measurements have been applied to characterize the reg- 
ularity of spatial distribution of two types of cones in 
agouti photoreceptor mosaics. The choice of the sum of 
angles instead of the hexagonality index was motivated 
by the fact that such systems involve relatively high levels 
of disorder. 

The agouti (Dasyprocta agout), is a hystricomorph 
rodent whose diurnal habit and well-developed visual 
streak yield an interesting model for comparative stud- 
ies of the visual system (e.g. 0,0). The organization 
of mosaics formed by short-wavelength-sensitive cones 
(S-cones) and middle- to long-wavelength-sensitive cones 
(M/L-cones) photoreceptors in the agouti's retina have 
been analyzed. For this purpose we used two polyclonal 
antibodies that have been shown to label S-cones (JH455) 
or M/L-cones (JH492) in a range of mammals, using im- 
munocytochemical methods largely according to the pro- 
cedures described elsewhere Q|. Using 40 x oil immer- 
sion objective, 25 fields taken with 250 x 250/im for M/L 
cones and 23 fields with 500 x 500/ito for S cones were 
acquired along the vertical-dorsal axis of the retina with 
the aid of an optical microscope (Eclipse E600, Nikon, 
Japan) equipped with a high-resolution video camera 
(Nikon 4500). The raw images of the cones were cap- 
tured at the level of the inner segments. For carrying 
out the analysis, x and y coordinates were identified by 
using the Scion Image Software (ScionCorp). 

Figure shows the average conformity ratios (ac-axis) 
and sum of angle differences (y-axis) obtained for each 
of the considered photoreceptor areas. It is clear from 
this figure that the two cone populations are clearly sep- 
arated by considering the sum of angle differences, while 
substantial overlap is observed for the conformity ratio. 
The obtained separation is even more definite than pre- 
vious results obtained for the same data set by using the 
lacunarity measurement of translational invariance 
The fact that the conformity ratio allowed the charac- 
terization of higher variability of the M cones suggests 
that additional insights can be obtained by using combi- 
nations of these two measurements. 

In order to investigate the effect of using successive 
shortest distances from each point, conformity ratios and 
coefficient of variations were obtained considering the 
second to the fifth shortest distances. Discriminant anal- 
ysis 0, El 0| was then applied in order to identify the 
contribution of each of such measurements to the overll 
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Measurement 


linearity 


sensitivity 


discriminative power 


Number of sides 
Coef. of variation 
Conformity ratio 
shortest distance 
sum of angle diffs. 
hexagonality index 


very smal 

high 

small 

high 

high 

medium 


very low 

almost constant 

higher for small perturbations 

almost constant 

almost constant 

slightly higher for small perturbations 


very low 

higher for small perturbations 
higher for large perturbations 
higher for large perturbations 
higher for small perturbations 
high and almost constant 



TABLE I: The linearity, sensitivity and discriminative power of the considered measurements of spatial order. 



separation between the two types of mosaics. The results 
indicated that the sum of angles is the only measurement 
among all possible linear combination of the considered 
features which allows all the mosaic types to be correctly 
identified. The use of conformity ratios and coefficients 
of variation for several shortest distances always implied 
in misclassifications of the mosaics. 



C. Characterization of Local Spatial Order 

Although we have so far concentrated on the charac- 
terization of the global properties of the spatial order in 
systems of points, another relevant problem involves the 
quantification of local order around each point. This is- 
sue is particularly important because, provided one can 
quantify the spatial order around a small neighborhood of 
each point, it becomes possible to study how such a mea- 
surement varies along the whole system of points. A par- 
ticularly interesting example of such cases, which often 
appears in complex systems, is the existence of multiple 
spatial domains characterized by varying levels of spatial 
order. This situation is illustrated in Figure [!3 including 
two simulated regions with different levels of spatial or- 
der. These domains, separated by an oblique border at 
the middle of the image, were obtained by perturbing a 
hexagonal lattice with 6 = 2 and 4, respectively. Such 
a relatively small variation of spatial order was inten- 
tionally imposed in order to impose a more demanding 
discrimination task on the measurements. 

Another interesting situation is characterized by the 
progressive variation, along space, of the spatial order. 
A typical example of such systems is the spatial distri- 
bution of photoreceptors in the mammals retina. Higher 
spatial order, actually almost perfectly hexagonal, is ob- 
served for those receptors near the central part of the 
retina, becoming less and less ordered as one moves to- 
ward the retina periphery. Figure ITU1 illustrates this type 
of point system, derived from a hexagonal lattice by 
adding perturbation values increasing linearly with the 
distance from the center of the system of points. 

Note that the conformity ratio and coefficient of vari- 
ation can not be used for quantifying local order, as the 



standard deviation of the shortest distances for a single 
point becomes equal to zero. Therefore, we limit the 
subsequent analysis to the hexagonality index, the sum 
of angle differences, and the first to third shortest dis- 
tances. Figure |H] shows the average ± standard deviation 
considering all individual points in each of the 50 realiza- 
tions of systems of points for each perturbation intensity. 
Except for the higher dispersions, these results are sim- 
ilar to those obtained by averaging globally over each 
realization shown in Figure [3J 

Figure |H1 presents the results of local order identifica- 
tion of the domains in Figure ??(a) by using the first (a), 
second (b) and third (c) shortest distances as well as the 
hexagonality index (d). In order to obtain the results in 
Figure^ the above measurements were obtained for each 
point and Bayesian decision theory (e.g. 0,0) was then 
applied in order to obtain the best threshold for sepa- 
rating between the two groups. The points identified as 
belonging to the higher order group are represented as the 
squares with holes in Figure^l Two particularly interest- 
ing results can be inferred from the obtained results: (i) 
the use of successive distances tend to reduce the correct 
identification of the regions and (ii) the hexagonality in- 
dex provided the best overall classification of points. The 
former effect suggests further insights about the fact that 
the consideration of further shortest distances does not 
tend to contribute significantly to the characterization of 
global order (see Section II V Bj l. 

Figure ^| presents the quantification of the local spa- 
tial order of the simulated retina by using the shortest 
distance (a) and sum of difference angles (b). The en- 
hanced sensitivity of the latter measurement is again cor- 
roborated by the more clearly defined gradient of hexag- 
onality values as one moves radially from the center of 
the structure towards its border. The higher uniformity 
of the spatial order of the central cells is also more clearly 
identifiable by using the sum of angle differences. 



V. CONCLUDING REMARKS 

The important problem of quantifying the spatial order 
of systems of points has been addressed. A series of re- 
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FIG. 4: Examples of simulated hexagonal lattices with perturbations 5 = 2 (a) and 4 (b). The hexagonality index is expressed 
by the gray levels at each Voronoi cell, with higher values indicated in light gray. Note the border elements in black, whose 
hexagonality is not considered. The minimum and maximum gray levels in each image are normalized between black and white 
for the sake of better visualization. 
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FIG. 5: The average ± standard deviation of the number of Voronoi cell sides (a), coefficient of variation of first shortest distance 
(b), conformity ratio of first shortest distance (c) first shortest distances (d), average of angle differences (e), and hexagonality 
indices (f), in terms of the perturbation intensity 5 considering all simulations. For the sake of better visualization, the standard 
deviations are shown to 5 times their original values. 
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quircmcnts expected from a good measurement of spatial 
order were identified, allowing the discussion and com- 
parison of a series of traditional (number of neighbors, 
nearest distances between pairs of points, conformity ra- 
tio and coefficient of variation of nearest distances) as 
well as a recently introduced (the hexagonality index) 
and a novel (the sum of angle differences) measurements. 

The potential of such measurements has been inves- 
tigated with respect to the quantification of global and 
local spatial order. In the former case, which involves 
assigning a single measurement to the whole system of 
points in order to characterize its order, the measure- 
ments were first compared with respect to simulated 
data, namely hexagonal lattices with progressive per- 
turbation intensities 5. The obtained results allowed 
the objective characterization of the linearity, sensitiv- 
ity and discriminative power of each considered mea- 
surement with respect to different levels of spatial or- 
der and hexagonality. The hexagonality index was veri- 
fied to account for good linearity and slightly enhanced 
sensitivity for higher spatial orders, with a nearly con- 
stant discriminative power. The sum of angle differences 
also yielded interesting properties including good linear- 
ity and discriminative power. The identification of the 
type of retinal mosaics in the agouti was also considered 
as a practical problem involving the estimation of spatial 
order. Because such systems involve relatively low spa- 
tial order, the sum of angle differences was used instead of 
the hexagonality index. Conformity ratios and coefficient 
of variations involving the first to fifth nearest distances 
were also calculated and used as subsidy for distinguish- 



ing between the two types of mosaics. Discriminant anal- 
ysis shown that the sum of angle differences was the only 
feature allowing perfect classification. 

The possibility to use the considered measurements 
to quantify local spatial order was also addressed with 
respect to progressively perturbed hexagonal lattices, 
yielding similar results as for the global analysis, but with 
substantially higher dispersions of measurement values 
for each perturbation intensity. Two important applica- 
tions of local order restimation, namely the identification 
of regions and analysis of systems of points involving gra- 
dients of spatial order, were also addressed. Such appli- 
cations further corroborated the good features of the sum 
of angle differences and hexagonality index as estimators 
of spatial organization. 

Such results imply that representative new information 
and insights about ordered biological systems can be ob- 
tained by revisiting previous investigations reported in 
the literature while considering a combination of sum of 
angle differences and hexagonality index measurements. 
Future works should include the extension of the angu- 
lar polygonality measurements to higher dimensions and 
applications to the characterization of texture and gene 
expression patterns. 
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FIG. 6: The histograms of number of Voronoi cell sides (first column), nearest neighbor distances (second column) and 
hexagonality indices (third column) obtained for 50 realizations of each of the configuration, with progressive perturbation 
intensities <5. 



10 



197p 

193- 

189- 

■ft 18S - 
t 

| 181 ; 

177- 
173- 
169- 
165- , 
161- 



♦ M 

x S 



x x 

X X 



. ♦♦♦♦ 

♦ 



157- 



0.21 0.25 0.29 0.33 0.37 0.41 0.45 0.49 

co^fficientafvariation 



[16] H. Waessle and H. J. Riemann. The mosaic of nerve cells 
in the mammalian retina. Proc. Royal Soc. B (London), 
200:441-461, 1978. 

[17] J. M. Ziman. Models of disorder. Cambridge University 
Press, 1979. 

[18] The property P may also be a function of a series of 
parameters s. The situation considered in this article as- 
sumes all but the single parameter s are fixed. 

[19] This type of statistical analysis involves the identifica- 
tion of the project from a high dimensional measurement 
space into a few dimensions so as to optimize the sepa- 
ration between the involved classes of objects. 



FIG. 7: Average conformity ratio and sum of angle differ- 
ences obtained for each of the cone mosaics. Note that the 
S and M/L cone populations can be completely separated 
by considering the sum of angle differences isolately, which is 
not possible in the case of the conformity ratio. On the other 
hand, the conformity ratio allowed the identification of a high 
variability of the M/L cone arrangements. 




FIG. 8: The average ± standard deviation of the number of Voronoi cell sides (a), first shortest distances (b), average of angle 
differences (c), and hexagonality indices (d), in terms of the perturbation intensity S considering all individual points in the 
simulations. 
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FIG. 9: Local order identification of a system of points involving two regions with distinct spatial order by using the first (a), 
second (b) and third (c) shortest distances and the hexagonality index (d). The squares with holes correspond to the points 
whose measurements were more likely to belong to the lower domain (higher spatial order.) 
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FIG. 10: Distribution of local order for the simulated retina by 
using the shortest distance (a) and the hexagonality index(b). 
The Voronoi cells belonging to the borders of the structure 
have been disconsidered. 



